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Abstract. In this paper, we prove existence and uniqueness of measure solutions for the Cauchy 
problem associated to the (vectorial) continuity equation with a non-local flow. We also give a 
stability result with respect to various parameters. 
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1. Introduction 
In this paper, we consider the system of nonlocal continuity equations 
f dtp' + div {p'V\t, X, r/**p)) = , t G M+ , x G 



(1.1) 



[p*(0)=p% iG{l,...,fc}, 

where the unknown p = {p^,...,p^) is a vector of measures, r/* = {rj^'^ , . . . ,7]^'^) is a vector of 
convolution kernels and we set ri'*p = {rj^'^ * p^, . . . ,if'^ * p^). For any time t ^ 0, if /xt G M.'^iW^) 
is a bounded measure on and r^t is a bounded function on W^, then the convolution is taken with 
respect to space only and is defined as usually as p-t * Vt = J^d Vti^ — v) ^Ptiu)- 

For example in [6], the authors consider the scalar conservation law where F is a nonlocal functional 
with respect to p. This equation stands for various models such as a sedimentation model, a supply- 
chain model, a pedestrian traffic model. For physical reason, in the following we are looking for 
positive solutions possibly with concentration, i.e. for any time t the solution has to be in Ai^ (M.'^)'' . 

Our goal here is to improve the results of [6] , not only by considering a system, but also lightening 
the hypotheses on V and tj. We prove here existence and uniqueness of weak measure solutions to 
(1.1). Let us introduce the following sets of hypotheses: 

(V): The vector field V{t, x, r) : M+ x M"^ x R'' — > ^d,k is uniformly bounded and it is Lipschitz 
in (x, r) G M"^ X M'^ uniformly in time: 

V G L°°(]R+ X M"^ X M'^') nL°°(R+,Lip(]R'^ x R'',Md,k)) ■ 

{rj): The convolution kernel r]{t,x) : R^ x M*^ ^ Myt is uniformly bounded and it is Lipschitz 
in X G M'^ uniformly in time: 

7] G L°°(M+ X M*^) nL°°(M+,Lip(M'^,Mfc)) • 
The main result of this paper is the following theorem. 
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1.1. Theorem. Let p € Ai^{W^)^. Let us assume that V satisfies (V) and rj satisfies (rj). Then 
there exists a unique solution p G L°°(M"*", A4+(M'^)'^) to (1.1) with initial condition p. 

We refer to Section 2 for precise notations and definitions, in particular for the notion of solution. 

1.2. Remark. Assume V satisfies (V) and r] satisfies (rj). Then Theorem 1.1 is completed by the 
following properties: 

• If p G Li(M'^, (M+)^) then p G <(fO(M+, Li(M'^, (M+ up to redefinition on a negligible set 
of times, and for all time t ^ 0, for all i G {1, . . . ,k} we have = II/^'IIli' 

• If p G (L^ n L~)(R'^, (]R+)'=) then p G (]R+, L°°(]R'', (M+)'')) and for all time t ^ 0, we 
have ||/o(i)llL°° ^ with C a constant dependent on Lip^(y), Lip,^(F), Lip^(r/) and 

WpWm- 

• Let p,(J ^ A4^(M'^)'^' such that for any i, \\p^\\j^ = ||^*||_a4- P ^^'^ solutions of 
(1.1) associated to the initial conditions p and a, then we have the estimate: 

n{pt,^t)^e'''Wiip,a), 

where K = Lip^.(F) + Lip^(y)Lip^.(?7)||p||_y^^ + Lip^(y)Lip^.(r?)||p||_y^ and Wi{pt,at) is the 

Wasserstein distance of order one between pt and at- 
These properties are described in Corollary 2.9 and in Proposition 4.2. The Wasserstein distance of 
order one is rigorously defined in Section 3. 

1.3. Remark. In Theorem 1.1 as well as in the other results of this papers, it is in fact sufficient to 
require that V^{t, x, r) is L°° in t, x and in r. Indeed, /9*ry* is uniformly bounded by ||/o||_y\^ II^IIl°° ~ 
M. Consequently, denoting Bm the closed ball of center and radius M in M*^, it is sufficient to 
have V' G L°°(M+ x M'^ x Bm)- 

Note also that, restricting the definition of V and r/ to the time interval [0, Tq], we obtain a solution 
defined on the same time interval. Consequently, we can as well ask only V and rj to be in time 
instead of L°°. 

The system (1.1) stands for a variety of models. Let us present first a macroscopic model of 
pedestrian traffic. In a macroscopic pedestrian crowd model, p is the density of the crowd at time 
t and position x and y is a vector field giving the speed of the pedestrian. According to the choice 
of V, various behaviors can be observed. Several authors already studied pedestrian traffic in two 
dimensions space (A^ = 2). Some of these models are local in p (see [3, 10, 16, 17, 19, 20]) ; other 
models use not only the local density p(t, x) but the entire distribution of p, typically they depend 
on the convolution product p{t) * rj (see [5, 6, 7, 8, 13, 21]) which represents the spatial average of 
the density. Within the framework of (1.1), we can study the models presented in [6, 7, 8]. In [6], 
the authors considered for V the expression 

V = v{p * r])v{x) , 

where f is a scalar function giving the speed of the pedestrians; r] is a convolution kernel averaging 
the density; and v{x) is a bounded vector field giving the direction the pedestrian located in x will 
follow. This model is more adapted to the case of panic in which pedestrians will not deviate from 
their trajectory and will adapt their velocity to the averaged density. Indeed, even if the density is 
maximal on a given trajectory, if the averaged density is not maximal, the pedestrians will push, 
trying all the same to reach their goal. This behavior can be associated with rush phenomena in 
which people can even die due to overcompression (e.g. on Jamarat Bridge in Saudi Arabia, see 
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[15]). A similar model was introduced in Piccoli &: Tosin [11, 21], where the authors instead of an 
isotropic convolution kernel, consider a nonlocal functional taking into account the direction in which 
the pedestrians are looking. 

In [6], the authors study the scalar case in the framework of Kruzkov entropy solutions. They ob- 
tained existence and uniqueness of weak entropy solutions under the hypotheses v € W^'°°(]R^", M"''), 
V € (W2'°° n W2'i)(M'^,M'^), and rj G {W'^'°° fl L1)(M^,]R). This result was slightly improved in [7] 
where, under the same set of hypotheses on v, v and r], the authors consider a system instead of a 
scalar equation and obtain global in time existence and uniqueness of entropy solutions. We recover 
these results with lighter hypotheses. Indeed, although we consider weak measure solutions, these in 
fact are unique and consequently coincide with entropy solutions when the initial condition is in L^. 

Another model of crowd dynamics that we recover consists in the coupling of a group of density 
p{t,x) with an isolated agent located in p{t). This can modelize for example the interaction between 
groups of preys of densities p and an isolated predator located in p. Such a model was introduced in 
[8] where the authors obtained existence and uniqueness of weak entropy solutions under very strong 
hypotheses. We recover here partially the results concerning the coupling PDE/ODE of [8]. Indeed, 
the measure framework allows us also to introduce particles/individuals through Dirac measures. 
For instance, let us assume that k = + ki such that /O^, . . . , p^^ are in fact functions belonging 
to L°°(M+,L^(M'^,R+)) and that J^i, . . .,5^k^ are Dirac measures located in p^{t), . . . ,p^^{t) G W^. 
Let us denote p = {p^,. . . , p^") G M^o and p = {p^, . . . .p^^) G ^dM- denote with (resp. 

rf) the vector fields (resp. kernels) associated to p, and with [/* (resp. A*) the vector fields (resp. 
kernels) associated to p. Note that 5p3 *X^'^{x) = X^'^ {x — p^). By definition of weak measure solution 
(see Definition 2.2), if pi G 'W^{[0,T],W^), the Dirac measures are satisfying, for any i G {1, . . . , ki} 

Piit) = [t, K(t), Pt*vliMt)), X'''{p^it)-plit)),...,y'''-{p,{t) -pkM 

which can be rewritten 

f{t) = <^' {t,p{t),pt*7^l{p\t))) . 

Consequently, in this case, the system (1.1) becomes 

V + div (^p'V'{t,x,r]i*pt, (A^'^(x-p?(t)))y) =0, iG{l,...,A;o}, 
p^{t) = (t, p{t), pt*rft{P^{t))) , j e{l,...,ki}. 

So we are coupling ODE with conservation laws. 

System (1.1) can also stands for models of aggregation, studied in [4] under weaker hypotheses 
admitting singular kernels. 

The system (1.1) comprised also a model of particles' sedimentation 

dtp + dx{{p * r]) p) = 

which has been introduced [22] and studied in [24] , where the author proved existence and uniquness 
of weak solutions with initial condition in L°°. 

Finally, a similar nonlocal model is the one the supply-chain model [1, 2], in which we consider 
the nonlocal term p{t, x) dx instead of a convolution product. This last model was studied for 
example in [9] with furthermore boundary conditions in x = and x = 1. 
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The proof of Theorem 1.1 is divided into two main steps. First, we prove some a priori properties 
of the solutions (see Section 2): mainly, we prove that the weak measure solutions of (1.1) coincide 
with the Lagrangian solutions of this system. Important consequences are the conservation of the 
regularity of the initial condition and the strong continuity in time in the case the solution is a 
function, as stated in Remark 1.2. 

Second, we prove the existence and uniqueness of Lagrangian solutions thanks to a fixed point 
argument (see Section 5). Indeed, introducing the set of probability measures endowed with the 
Wasserstein distance of order one, we are able to prove a stability estimate with respect to the 
nonlocal term (see Section 4). The technique used there is quite similar to the one of Loeper [18], 
who studied the Vlasov-Poisson equation and the Euler equation in vorticity formulation. 

This article is organized as follows: in Section 2 we define the two different notions of solution 
and prove that they coincide. In Section 3 we give some useful tools on optimal tranport; in Section 
4 we prove an important lemma giving a stability estimate and in Section 5, we give the proof of 
Theorem 1.1. 

2. Notion of solutions 

2.1. General notations. Let d G N be the space dimension and G N be the size of the system. In 
the following, M^^^ is the set of matrices of size dx k with real values and is the set of matrices 
of size k X k with real values. We denote by A^(R'^) (resp. the set of bounded (resp. 
bounded and positive) measure on M.'^ and by 'P(M'^) the set of probability measures on M'^, that is 
the set of bounded positive measures with total mass 1. 

In the following the Lipschitz norms with respect to x or r are taken uniformly with respect to 
the other variables. That is to say, for example: 

Lip,(y)= sup {Up,{V{t,-,r))}. 

Let us also underline that in the computations, we considered the norm 1 on the vectors in 
M'^. When considering another norm, a constant depending on k appears in the various estimates. 
Similarly, if p = {p^, . . . , p^) G ^A{Mf^)^ we define the total measure of p as \\p\\j^ = \\p^\\j^ + . . . + 

ll^'^II.M- 

The space L°°([0, T], 7W"'"(M'^)) consists of the parametrized measures p = {pt)t£[o,T] such that, for 
any (/) G '^^'^(M'^, M), the application t J^a <ipt{x) is measurable and such that ess supigp ^] ll^t|lA4 < 
oo. 

2.2. Weak measure solutions. We say that p G Ij°°{[0,T],A4^(R'^)'') is a weak measure solution 
of (1.1) with initial condition p G A4~^{W^)'' if, for any i G {1,...,A;} and for any test-function 

(/) G '^c°°(] - oo, r[xM'^, R) we have 

/ / [dt^ + V\t,x,p*r]')-\7(l)]dpl{x)dt+ [ (l){0,x)dp\x)=0. 
Jo JW JRd 

2.3. Remark. A priori for weak measure solutions of the continuity equation dtp + div (pb) = 0, with 
a given vector field b, we have only continuity in time for the weak topology (see [12]), that is to say, 
for all i G {1, . . . , k}, for all (p £ (^'^i ^)' the application t i— )• (t>{x) dp\{x) is continuous, up to 
redefinition of pt on a negligible set of times. 
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In the case of the system (1.1), we have a gain of regularity in time when the initial condition is 
a function in L'^{R'^, (R+)^) (see Corollary 2.9). 

2.4. Push-forward and change of variable. When ^ is a measure on Q and T : il. ^ il.' a 
measurable map, we denote Tj/i the push-forward of fi, that is the measure on Q' such that, for every 
(/) G ^0(J1',M), 



Jn' Jn 



If we assume that fi and v = T^fi are absolutely continuous with respect to the Lebesgue measure 
so that there exist f,g e such that dfi{x) = f{x)dx and du {y) = g{y)dy, and that T is a 
Lip-diffeomorphism, then we have the change of variable formula 

f{x) = g{T{x))\det{VT{x))\. (2.1) 

Besides, we denote P^;. : M'^ x M*^ — )• M'^ the projection on the first coordinate; that is, for any 
(n, v) G M'^ X M'^, Fx{u, v) = u. In a similar way, Pj^ : M'^ x M*^ — > M"^ is the projection on the second 
coordinate; that is, for any {u,v) G M"^ x W^, Wy{u,v) = v. 



2.5. Lagrangian solutions. We say that p G L°°([0, T], A^+(M'^)'^) is a Lagrangian solution of 
(1.1) with initial condition p G 7W~^(M'^)'^ if, for any i G {1,...,A;}, there exists an ODE flow 
X* : [0, T] X M'' ^ R'^, that is a solution of 

^(t,x) = {t,X\t,x),pt*4{XKt,x))) , 
X\Q,x)= x; 

and such that p\ = Xl^p' where XI : W'- ^ is the map deflned as Xl{x) = X^{t,x) for any 
{t, x) G [0, T] X R"^. 

2.6. Remark. Assume V satisfies (V) and t] satisfies (r/). Then, for any p G L°°([0, T], 7W+(M'^)''), 
the vector fields b = V{t,x, pt*r]t) Lipschitz in x and 

LiPx(^) ^ Lip^(^) + ^^Pr{y)^^Px{v)\\pt\\M ■ 

Consequently, if ||/3t||;\4 is uniformly bounded, the ODE flow X^ above is always well-deflned, for a 
fixed p. 

If p G L^(]R'^, IR+), then the push-forward formula (2.1) becomes, for a.e. {t,x) G M+ x W^, 

p\t, X\t, x)) = p\x) exp {- div V' (r, X^t, x),pr* Vr{X'{T, x))) dr^ . (2.2) 

We now show that the two notions of solution in fact coincide. 

2.7. Theorem. If p is a Lagrangian solution of (1.1), then p is also a weak measure solution of 
(1.1). Conversely, if p is a weak measure solution of (1.1), then p is also a Lagrangian solution of 
(1.1). 
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Proof. 1. Let p be a Lagrangian solution of (1.1). Let us denote 6* = V^{t,x, p*r]^) and let X* be 
the ODE flow associated to b\ Then, for any (/> G '^^{] - oo, r[xM'^, M), we have 

{dtcpit, x) + b\t, x) ■ V(p{t, x)) dptix) dt 



{dt<P{t,Xl{x)) + b\t,Xi{x)) ■ Vcl>{t,Xl{x))) dp{x)dt 
dp(x)dt 

(f>{T,X\T,x))dp{x)- [ (f>{0,x)dp{x) = - [ (p{0,x)dp{x 



which proves that p is also a weak measure solution. 

2. Let p be a weak measure solution of (1.1). For any i € {1,...,A;}, let us denote ¥{t,x) = 
V^{t,x, p*r]^). Let a be the Lagrangian solution of the equation 

dta' + div {a%') = , a'{0,-)=p\ (2.3) 

which exists and is unique since 6* is Lipschitz as noted in Remark 2.6. Then, arguing similarly as 
in point 1, a is also a weak measure solution to (2.3). Denoting = — o"*, we obtain that ti* 
is a weak measure solution of the equation dfu"^ + div (n* ¥) = with initial condition ti*(0, •) = 0. 
Consequently, for any (/) £ 'tf^(\ - oo,T[xR'^,R), 

[ [ {dt^ + b'{t,x)-V(l))dutdt = 0. 
Jo 

Let V € "tf^i] - oo, r[xM'^, M). Since 6* G L°°([0, T] x R'^, R'^) is Lipschitz in x, by computation along 
the characteristics, we can find (/) G '^^Q - oo,r[xM'^,M) so that ip = dt(j) + 6*(t, x) • V(j). Hence, for 
any ip G ^c(\ ~ TfxM'', M), we have Jj^a ^ dut dt = 0, which implies n = a.e. and p = a a.e. 
Consequently, we have also ¥{t,x) = V^{t,x,a*if ), and cr = p is finally a Lagrangian solution of 
(1.1). □ 

2.8. Definition. As a consequence of the previous theorem, in the following we simply call solution 
of (1.1) a weak measure solution or a Lagrangian solution of (1.1), that in fact coincide. 

It is now possible to prove some of the properties given in Remark 1.2. 

2.9. Corollary. Assume that V satisfies (V) and ij satisfies (rj). Let p G L°°([0, T], A1+(M'^)'') be a 
solution to (1.1) with initial condition p G M.'^iR'^)^ . 

• Ifpe L^iR"^, {R+)''). Then we have p G 'r°([0, T], L1(M'^, (1R+)*^)) and for all time t G [0,r], 
alii G {1,...,A;}, UpHOUli = II^IIli- 

• If furthermore p G (L^ nL°°)(]R'^, (R+)''), then for allte [0, T] we have p{t) G L°°(IR'^, {R+f) 
and we have the estimate 

llp(t)llL- ^ llpllL-e^*, 
where C depends on \\p\\j^, V and r}. 

Proof. Let phe a solution of (1.1) with initial condition p G L"'^(IR'^, (M'^)'^). According to Definition 
2.8, p is a Lagrangian solution associated to a fiow X and we have immediately that ||p||li = ||p(i)llLi- 



SYSTEM OF CONTINUITY EQUATIONS WITH NON-LOCAL FLOW 



7 



Besides, as U'{t,x) = V'^{t,x, p*if) € L°°([0, T] x M'^,R'^) is bounded in t and Lipschitz in x, then 
Xi G Lip(]R'^,M'') and we can use the change of variable formula (2.2). If p G L°°(]R'', M^), with 

C = Lipjy) + Lip,,(y)LipJr/)||p|t^ , 

we obtain the desired L°° bound and p{t) G L°° for all t G [0,T]. 

The continuity in time can be proved directly by estimating ||pt — PsIIli using Egorov Theorem. 
This computation is straightforward although a bit long so we prefer to omit the details. 

Besides, note that the continuity in time is also ensured by the results of DiPerna & Lions [14, 
Section 2. II] and the notion of renormalized solutions. 

□ 



3. Some tools from optimal mass transportation 
Let us remind the definition of the Wasserstein distance of order 1. 

3.1. Definition. Let p, v be two Borel probability measures on W^. We denote H (^, z^) the set of 
plans, that is the set of probability measures 7 G A4^{W^ x W^) such that IPxjjT = P ^-nd = v- 
We define the Wasserstein distance of order one between p and u by 



Wi{p,v) = inf / \x - y\(i-i{x,y) . (3.1) 

Let p = {p^, . . . , p^), a = (cj^, . . . , fj'^) be two vectors such that p^, . . . ,p^ and a^, . . . are Borel 
probability measures on M*^. We define the Wasserstein distance of order one between p and a, 
denoted Wi{p,a), as 

k 

Wi{p,a) = Y,Wi{p\a'). (3.2) 

i=l 

3.2. Remark. By [23, Theorem 1.3], for any p,i' £ 'P(M'^), there exist a plan 70 G H(/u, z/) realizing 
the minimum in the Wasserstein distance so that 



Wi{p,i')= \x - y\d-fo{x,y) . 

3.3. Remark. Let p G A4~^(M.'^)^ be a probability measure ; and let X, y : M'^ ^ M'^ be mappings 
such that f = X^dp and g = dp. Then, the probability measure 7 = {X, Y)^ dp satisfies Fxj7 = /, 
Pyjj7 = g and so 

Wi{f,g)^ [ \x-y\dj{x,y)= [ \X-Y\dp{x). 

3.4. Proposition (cf. Villani [23, p. 207]). Let p^v he two probability measures. The Wasserstein 
distance of order one between p and v satisfies 



Wi{p,v) = sup / (l){x) {dp{x) — dv{x)) . 

Lip{0)^l Jm.<i 
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4. The main stability estimate 

In the following we consider probability measures instead of bounded positive measures. This is 
not a real restriction since we pass from one case to the other just by a rescaling. 
Before giving a stability estimate in Proposition 4.2, we prove a technical lemma. 

4.1. Lemma. Let V satisfy (V) and rj satisfy (r]). Let r, s € V{M.'^)^. For any i G {1, . . . we 

have the following estimate 

\\V\t,x,r*7^l)-V\t,x,s*iil)\\^^ ^ Lip,(nLip,(r?^)#i(r,s). 

In the previous lemma, the quantity #i(r, s) could be infinite. If we restrict ourselves to bounded 
positive measures with first moment finite, then the quantity above is always finite. 

Proof. The proof follows from Proposition 3.4 on the Wasserstein distance. Note first that in the 
case Lip(r/*'-') = then rj'''^ is constant and we have (r-' — s^) *if'^{x) = ^ Lip{r]^'^)Wi{r^ , s^). Now, 
in the case Lip(?7*'-') ^ 0, thanks to Proposition 3.4, we have 



{r^ - s^) *if'\x) = I Tf'^{x-y){dr^{y)-(ls\y)) 

= Uv{lH)j 'p^^{dr^{y)-ds^{y)) 



^Lip(?f'^) sup / (l){y){dr^y) -ds^iy)) =Uj){if'^)Wi{r^,s^] 

Lip((/))^l JRd 

As we obtain the same estimate for —{r^ ~ s-^) * Tf'^{x), we can conclude that 

\\V\t,x,r*rt)-V\t,x,s*rt)\\^^ ^ Uv,{V%\{r - s)*r^W^^ 

^ UVr{V')Uv{rt)Wi{r,s). 



□ 



Let r,s G L~([0, T], P(1R'^)^). We want to compare the following equations, in which the nonlocal 
has been replaced by fixed applications, so that the system is made of decoupled equations. 

for ah i G dtp' + diw {p'V'{t,x,ri'*r)) = {) , p\Q^.)=p\ 

\ . . . \ (4.1) 

for alH G {l,...,/c} 9(0-* + div (cj*[/*(t, x, z^**s)) = , a\^,-)=a\ 

4.2. Proposition. Assume V, U satisfy (V) andrj, v satisfy ("q). Let p, a be two probability measures 
such that for any i, = ||^*||_vc r,s ^ L°°([0, T], 'P(]R'^)^). // p and a are Lagrangian 

solutions of (4-1) associated to the initial conditions p and a, then we have the estimate: 



HpuTt) ^ e^''Wi{p,a) + Cte^' 



sup #L(n, St) + \\ri - i^IIloo + \\V - U\\^o 
fe[o,T] 



(4.2) 



where C is a constant depending on Lip^.(y), Lip^(y), Lip^.(r/) and \\p\\_x4- 
Furthermore, in the special case r = p and s = a, we get: 

n{puat) ^ e'''n{p,^) + Kte'''[\\r^-u\\^^ + \\V -U\\^^] , (4.3) 

where Kis a constant depending on Lip^(F), Lipj,(y), Lip^(?7) and \\p\\j^- 

Note that the estimate above comprises the case y^(/0, a) = oo. 
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Proof. Let p, a be two Lagrangian solutions to the Cauchy problem for (1.1) with initial conditions 
p and a respectively. Let X, Y be the associated ODE flows. For any t E [0, T], we define the map 
Xi M y/ : M*^ X R'^ ^ M'^ X by 



Xi^Y^ix,y) = {Xtix),Y;iy)) 



for any (x, y) € M x 



Let 7q € E{p^,a^) so that Px-^Tq = and ^y^% = cr*. Let us define the probability measure 
7| = {XI N Yi)i%. Then, F^^f, = pi and Fy^t = al so that ji G S 
Let i? > 0, we define, for t ^ 



Quit) = E / 



_^ JXl[Ba)xY,\Ba) 
Note first that is Lipschitz. Indeed, let t, s ^ 0, then we have 



{y)\ d7^(x,y) . 



\QR(t)-QR{s)\ ^ 



/ (|x,^(x)-y/(y)|-|xi(x)-y;(y)|)d7^(x,y) 
E / |^l(x)-y/(y)-x:(x) + y;(2/)|d7j(x,2/) 

E / (|Xl(x)-X:(x)| + |y/(y)-y;(y)|)d7j(x,y) 

4=1 JBrXBr 



i=l JBrxBh 

^ (II^IIl- + llf^IlL-) ll7o|lAi|i-«l- 



(4.4) 



Let us assume that #i(/3, o") < oo, otherwise the thesis is trivial. Then, by Remark 3.2, for all 
i € {1, A:}, we can find a bounded positive measure 7q G H {ff, ct*) so that 



Wiip\a' 



Consequently we have, for any 0, Qr{0) ^ Wi{p,a). Hence, using (4.4), for any t ^ 0, we have 

QB.{t) < Qi?(0) + (||^7||L^ + ||y||Loo)||70|lA4* 

5^ Wiip,a) + i\\U\\j^^ + ||1/||loo)||7o|Ia4*- 

Thus, for any t ^ 0, Quit) remains finite when i? — ?> oo and since R i— )■ Qnit) is increasing with 
respect to R, we can define Q{t) = lim/j_!,oo Qnit)- 
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Let US now consider Q. The same computation as in (4.4) ensures that Q is Lipschitz so we can 
differentiate for almost every t and obtain 

k 



i=l -'K'^xIR'' 

^ E / f|^^(t,^l(x),r,*r?j(Xl(x)))-r(t,y/(y),r,*r/j(X,^(x)))| 

+|r {t,Yiiy),n*4mx))) - V' (t,y/(2/),ri*r?i(y/(y)))| 

+ |r {t,Y:{yln*4{Yl{y)))-V' (t, y/(y), s<*r?j(r/(y))) | 



Note that 

\V\t, y, n * - V\t, y, n * vHv))] ^ Lip,(y^)Lip(ri*r?j)|x - y\ 

^ Up,iV')\\rt\\j^Up{rj')\x-y\. 

Using Lemma 4.1 we obtain 

Q'it) ^ E/ (Lip,(n+Lip,(ni|ri||^Lip(r/^))|x;(x)-y/(y)|d7^(x,y) 

+ I UMV')Uv{v')n{ruSt)d^i{x,y) 

+ / Up,.{V')\sMvl-i^l)\dfo{x,y)+Y^ / ||r-C/*||Lood7^(x,y) 

i=l •^K'^xR'* ^^-^ JK'^xR'* 

^ C[Q(t) + #i(n,st) + ||r/-i/||Loc + ||?7-y||Loo] . (4.5) 
Taking the sup in time of Wi{rt, st) on the right-hand side and applying Gronwall Lemma, we get 

Q{t) ^ e^*Q(0) + Cte^* ( sup #i(r,, s^) + \\v - i^IIl- + 11^ " ^^IIl- I • 

\re[o,t] J 

Note now that, thanks to remark (3.3), for any t ^ 

mpt,at) ^ Q{t). (4.6) 

Furthermore, we have chosen 70 in an optimal way thanks to Remark 3.2 so that Q{0) = Wi{p,a). 
Hence we obtain, for any t € [0, T]: 

nipt, at) ^ e^*#i(p, a) + Cte^* ( sup #i(r,, Sr) + Wv - + \\U - V\\j^^ ] , 

\re[o,t] J 

which is the expected result (4.2). 
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In the particular case r = p and s = o", applying (4.6) to (4.5) we obtain 

Q'{t) ^ 2CQ{t) + C{\\i]- z^IIloo + \\U - FIIloo) . 

Applying Gronwall Lemma, we finally obtain Q{t) ^ e^'^*(5(0) + Cte'^'^^ {\\r] — v\\i^oo + \\U — V\\^oo), 
which is (4.3). 

□ 

5. Proof of the main theorem 

The proof of Theorem 1.1 is based on the following idea: let us fix the nonlocal term and, instead 
of (1.1), we study the Cauchy problem 

dtp + dw{pV{t,x,r*7])) = 0, p{0)=p, (5.1) 

where r is a given application. We consider here probability measures. In the more general case 
of positive measures with the same total mass, by rescaling we are back to the case of probability 
measures. 



Let us introduce the application 



=S : ^ ^ ^ ^ ^ , (5-2) 



where we consider the space = Ij°°{[0,T],V{W^)'') for T chosen in such a way that: 

(a) : The space ^ is equipped with a distance d that makes ^ complete: for p^u ^ we 
define 

d{p,v) = sup Wi{pt,ut) ■ 
te[o,T] 

(b) : The application ^ is well-defined: the Lagrangian solution p E to (5.1) exists and is 
unique (for a fixed r). Indeed, let Xt be the ODE flow associated to V{t,x,rt*r]t), then we 
can define pt = Xf^p. Since p is a positive measure, then so is pt- 

(c) : The application ^ is a contraction: this is given by Proposition 4.2. Indeed, let r,s in 
L°°{[0,T],M+{R'^)'') and denote p = ^(r), a = =2(s) the associated solutions to (5.1). Note 
that p and a have the same initial condition. Thanks to Proposition 4.2, we obtain the 
contraction estimates 

sup Wiipt, at) ^ CTe^^ sup Wiin, st) , 

[0,T] [0,T] 

where C depends only on Lip^{V), Lip^(y), Lip^{r]) and 

Hence, for T small enough, by the Banach fixed point Theorem we obtain existence and uniqueness 
in ^ of a Lagrangian solution to (1.1) for t G [0, T]. As ||/5t||_a4 ~ IIpIIai coefficient C does 
not depend on time and we can iterate the procedure. Thus we have existence and uniqueness on 
[0,+oo[. 

Observe that uniqueness can be also obtained directly by the stability estimate (4.3) in the par- 
ticular case = , rf = , p = a . 
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